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HIGHLIGHTS 


•  We  simulate  electrolyte  transport  dynamics  in  the  compressed  electrode  of  LIB. 

•  The  anode  has  lower  wettability  than  the  cathode. 

•  The  particle  shape  has  a  strong  effect  on  the  wettability. 
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Using  the  multiphase  lattice  Boltzmann  method  (LBM),  the  electrolyte  transport  dynamics  in  the  two- 
dimensional  electrode  structure  of  a  lithium-ion  battery  are  simulated.  The  effect  of  the  compression 
ratio  of  a  porous  electrode  on  wettability  is  explored  with  respect  to  variations  of  porosity  and  particle 
shape.  The  electrolyte  distribution  in  the  electrode  and  the  electrolyte  saturation  profile  are  examined  in 
order  to  evaluate  the  wetting  capability  at  various  compression  ratios.  The  results  show  that  wettability 
in  the  electrode  decreases  as  the  compression  ratio  increases.  In  a  highly  compressed  electrode,  the 
through-plane  permeation  of  liquid  electrolyte  is  small.  Thus,  the  electrolytes  are  mainly  observed  at  the 
interface  of  the  electrode  and  separator.  The  anode  has  lower  wettability  than  the  cathode  due  to  the 
deformation  of  particle  shape  during  the  manufacturing  process.  Therefore,  particle  shape  has  a  strong 
effect  on  wettability.  The  two-dimensional  LBM  approach  used  in  this  study  characterizes  the  electrolyte 
transport  phenomena  inside  the  electrode  and  allows  us  to  compare  the  wettability  between  the  cathode 
and  anode  at  various  compression  ratios. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  wetting  capability  of  an  electrolyte  is  a  critical  issue  for  the 
commercialization  of  lithium-ion  batteries  (LIBs)  in  electric  vehicle 
(EV)  and  hybrid  electric  vehicle  (HEV)  applications.  The  electrode 
needs  to  be  fully  wetted  with  the  electrolyte  to  maintain  high 
utilization  of  electrode  capacity.  Insufficient  electrolyte  wetting  in 
the  electrode  causes  irregular  reactions  of  the  electrochemistry, 
and  deteriorates  LIB  performance.  This  phenomenon  causes 
extrusion  of  lithium  metal  that  may  result  in  a  serious  problem  in 
relation  to  the  safety  issue.  Moreover,  inferior  wetting  with  the 
electrolyte  accelerates  the  degradation  of  LIBs,  resulting  in  a 
shortened  cycle  life. 
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Recent  research  efforts  for  EV  and  HEV  applications  have 
considered  high  energy  and  high  power  densities.  According  to  the 
United  States  Advanced  Battery  Consortium  (USABC),  long  term 
goals  for  energy  and  power  densities  are  300  Wh/L  and  600  W/L, 
respectively  [1  .  To  achieve  these  goals,  a  high  packing  density  of 
the  electrode  is  necessary.  On  the  other  hand,  this  may  lead  to 
insufficient  electrolyte  wetting  in  the  electrode  by  lowering  the 
porosity  and  increasing  pore  blockage.  Yu  et  al.  [2]  suggested  that 
thick  and  dense  electrodes  degrade  discharge  capacity  rapidly  at 
high  current  rates.  This  is  believed  to  be  due  to  the  insufficient 
wetting  in  the  electrode  caused  by  high  packing  density.  The  main 
objective  of  the  present  study  is  to  investigate  the  effect  of  elec¬ 
trode  compression  on  the  wettability  in  the  electrodes. 

A  number  of  experiments  have  been  conducted  regarding 
wettability  in  batteries.  Wu  et  al.  [3]  investigated  the  wettability  in 
the  electrodes  of  an  LIB.  Their  experimental  results  showed  that  an 
increase  in  surface  tension  reduces  electrolyte  spreading  in  the 
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pores.  They  reported  that  the  wetting  in  the  electrode  depends  on 
electrolyte  penetration  and  spreading  which  are  determined  by 
viscosity  and  surface  tension.  Zguris  [4]  and  Culpin  5]  examined 
the  wettability  of  the  separator  for  valve-regulated  lead/acid 
(VRLA)  batteries.  They  suggested  that  the  fluid  absorption  proper¬ 
ties  of  the  separator  are  an  important  battery  design  consideration. 
Manev  et  al.  [6]  examined  that  the  effect  of  preliminary  treatment 
and  technique  on  the  charge/discharge  behavior  of  graphite.  They 
reported  that  a  higher  compacting  pressure  on  the  electrode  leads 
to  a  decrease  in  the  specific  capacity.  Chu  et  al.  [7  investigated  that 
the  influence  of  compaction  on  the  porosity  and  electrochemical 
properties  of  positive  electrode.  They  suggested  that  the  wettability 
is  predominant  at  high  C  rate  along  with  densification.  Although 
many  efforts  have  been  made  regarding  wettability,  no  in-situ 
visualization  method  has  been  reported  that  can  detect  the  elec¬ 
trolyte  distribution  inside  the  electrode.  More  recently,  Lee  et  al.  [8] 
identified  the  wetting  mechanism  in  the  porous  electrode  of  an  LIB 
using  the  two-dimensional  lattice  Boltzmann  method  (LBM).  They 
demonstrated  the  dynamic  behavior  of  the  liquid  electrolyte  in  a 
porous  electrode  as  well  as  wetting  phenomena  in  a  microstruc¬ 
ture.  Their  results  demonstrated  that  the  LBM  can  be  effectively 
used  to  estimate  the  liquid  electrolyte  behavior  in  a  porous  elec¬ 
trode  by  giving  a  realistic  pore-scale  dynamic  simulation. 

The  LBM  has  been  widely  used  to  study  the  viscous  coupling 
effects  in  immiscible  two-phase  flow  in  porous  media.  Unlike  the 
traditional  computational  fluid  dynamics  (CFD)  method  which 
solves  the  Navier-Stokes  equation  at  a  macroscopic  scale,  the  LBM 
simulates  fluid  flows  by  tracking  the  evolution  of  a  single-particle 
distribution.  The  LBM  is  computationally  accurate  and  efficient 
method.  Due  to  this  characteristic,  the  LBM  can  separate  the 
immiscible  fluids  without  tracking  interfaces  while  sharp  interfaces 
are  maintained  automatically.  Therefore,  the  LBM  can  be  applied 
successfully  to  simulate  immiscible  flows  such  as  bubble  collision 
and  rising,  and  wetting  and  spreading  phenomenon.  The  LBM  is  a 
parallel  method  because  all  the  information  transfer  is  local  in  time 
and  space  [9,10].  Using  modern  computers  with  massively  parallel 
processors,  the  LBM  has  become  a  promising  computational 
method  for  studying  various  complex  systems  11]. 

This  paper  develops  a  numerical  simulation  using  the  LBM  to 
investigate  the  effects  of  the  compression  ratio  (CR)  on  wettability. 
The  two-phase  intermolecular  potential  model  proposed  by  Shan 
and  Chen  [9,10]  was  used  to  examine  the  microscopic  behavior 
occurring  at  the  pore  scale  in  a  complex  porous  structure.  Based  on 
our  previous  work  [8]  in  which  the  wetting  mechanism  in  a  porous 
electrode  was  investigated,  the  present  study  extends  the  wetting 
phenomena  in  both  the  cathode  and  anode  by  considering  various 
compression  ratios.  The  changes  of  particle  shape  are  taken  into 
consideration  in  the  anode,  because  the  anode  material  can  un¬ 
dergo  plastic  deformation  during  the  manufacturing  process. 
Cross-sectional  liquid  electrolyte  distribution  and  saturation  pro¬ 
files  are  presented.  This  study  focuses  on  a  comparison  of  wetting 
capability  between  the  cathode  and  anode. 


2.  Model  description 

2.1.  LBM  model 

The  LBM  originated  from  lattice  gas  automata  (LGA)  and  is 
numerically  derived  from  the  continuous  Boltzmann  equation 
discretized  in  time  and  phase  space  [12  .  In  the  LBM,  the  particle 
population  in  the  LGA  model  is  replaced  by  the  distribution  func¬ 
tion.  The  time  variation  of  the  distribution  functions  is  calculated 
by  performing  simple  law  of  collision  and  streaming  steps.  In  this 
study,  we  implemented  the  lattice  Bhatnagar-Gross-Krook  (LBGK) 


model  for  the  collision  term  approximation.  The  discrete  particle 
distribution  functions  fi  are  written  as 


fi(x  +  e,At,  t  +  At) 


At 


fi(x,  t)  -f,eq(x,  t) 


+  Sj(x,  t) 


where  is  the  probability  of  finding  a  particle  in  the  ith  velocity 
e*  at  (x,t),  and  At  is  the  time  step,  t  is  the  dimensionless  relaxation 
time  that  is  related  to  the  kinetic  viscosity  v  and  controls  the  ten¬ 
dency  of  the  system  to  relax  the  local  equilibrium.  S/(x,t)  is  a  source 
term  added  to  the  standard  lattice  Boltzmann  equation.  Various 
fluid  dynamics  phenomena  can  be  expressed  by  using  different 
forms  of  the  source  terms.  /Jeq  is  the  equilibrium  distribution 
function  which  is  a  specially  discretized  Maxwell-Boltzmann 
function  with  a  Taylor  expansion  12]  and  is  given  as  follows: 


pWf 


( ueq)2\ 

2  c\  J 


where  cs  is  the  speed  of  sound  in  a  liquid  electrolyte  which  is  c/3;  c 
denotes  the  ratio  of  lattice  space  to  the  time  step,  which  is  c  =  Ax/ 
At;  and  w/  is  the  weight  factor,  which  is  4/9  for  i  =  1,  1/9  for 
i  =  2, 3, 4, 5  and  1/36  for  i  =  6, 7, 8, 9.  In  this  study,  two-dimensional 
nine  velocities  (D2Q9)  model  was  used  on  the  two-dimensional 
square  lattice  as  shown  in  Fig.  1.  By  applying  the  D2Q9  model,  the 
discrete  velocity  e\  can  be  written  as 


[e\  ,  e2>  e3>  e4>  e5>  e6>  e7>  e8>  e9] 


0  10-1  0  1-1-1 
0  0  1  0-11  1  -1 


In  a  Newtonian  fluid,  the  fluid  density  p  and  equilibrium  velocity 
ueq  can  be  calculated  using  the  following  equations: 


9 


p  = 

i=  1 


Ueq  =  UH -  (5) 

p 

where  u  is  a  local  velocity  given  by  u  =  ^afa^a/p*  F  is  the  total  force 
acting  on  the  fluid  that  is  comprised  of  inter-particle  force  Tint> 
adhesion  force  Fa ds  between  a  liquid  and  gas  with  a  solid  phase,  and 
external  force  Text*  In  this  study,  it  is  assumed  that  Tex  t  has  a  uni¬ 
form  value  in  a  steady  body.  Fmt  and  Facjs  were  calculated  using  the 
following  equations: 


Fig.  1.  Lattice  structure  of  two-dimensional  nine  velocity  (D2Q9)  model. 
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Fig.  2.  Schematic  of  phase  interfaces. 


Fig.  4.  Reconstruction  of  cathode  electrodes  when  the  compression  ratios  are  (a)  0%, 
(b)  10%,  (c)  20%,  (d)  30%  and  (e)  40%. 


9 

Fint(*.  t)  =  -Gi(x,  f)  J2  Wji(x  +  CjA t,  t)c,  (6) 

l'=l 

9 

Fads (x,t)  =  -G\p(p(x,  t))  ^2  Wj\p(p(xw)s(x  +  e,At,  t)e,-  (7) 

l'=l 

where  G  is  Green’s  function,  which  controls  the  strength  of  the 
interparticle-force,  and  the  signs  of  G  determines  whether  it  is 
attractive  or  repulsive  [8  .  f  is  a  function  of  density  in  the  following 
empirical  form:  \j/(p)  =  po[l  -  exp(-p/po)],  where  po  is  constant,  s  is 
an  indicator  function:  1  for  a  solid  domain  node  and  0  for  a  fluid 
domain  node.  p(xw)  is  the  density  of  solid  phase,  and  provides  a  free 
parameter  with  which  to  tune  different  wall  properties  indepen¬ 
dently  from  the  real  density  value.  The  desired  contact  angle  can  be 
derived  by  controlling  the  value  of  p(xw)  13].  This  allows  the 
different  contact  angles  for  various  wall  conditions  to  be  easily 
imported. 

This  proposed  model  adopts  the  Shan-Chen  pseudo-potential 
model  which  is  one  of  the  most  popular  multiphase  models.  The 
pseudo-potential  model  allows  the  separation  of  two  immiscible 
fluids  if  the  interaction  force  between  the  fluids  is  larger  than  a 
threshold  value.  The  multiphase  interface  includes  fluid-fluid 
interface  and  fluid-solid  interface  as  shown  in  Fig.  2.  A  fluid-fluid 
interface  is  formed  between  the  two  immiscible  fluids  and  fluid- 
solid  interface  is  formed  between  fluid  and  solid.  According  to 


(b) 


Fig.  3.  Schematic  of  roll  pressing  process  for  (a)  cathode  and  (b)  anode  electrodes.  The 
corresponding  images  of  the  particle  are  presented.  Note  that  the  cathode  particles 
maintain  spherical  shape,  whereas  the  anode  particles  are  deformed  to  elliptic  shape 
after  roll  pressing  process. 


Fig.  5.  Reconstruction  of  anode  electrodes  when  the  compression  ratios  are  (a)  0%,  (b) 
10%,  (c)  20%,  (d)  30%  and  (e)  40%.  The  deformation  of  particles  is  considered. 

Srivastava  et  al.  14],  the  Shan-Chen  multiphase  model  and  the 
multi-component  LBM  model  shows  good  agreement  with  the 
sharp  interface  model. 

2.2.  Reconstruction  of  porous  electrodes 

The  commercially  available  cathode  materials  are  LiCo02  (LCO), 
LiMn204  (LMO),  Li(NiCoMn)02  (NCM)  and  LiFeP04  (LFP).  These 
materials  are  ceramics,  which  are  hard  and  strong  in  compression. 


Fig.  6.  Predicted  porosity  at  various  compression  ratios  with  error  bars  which  repre¬ 
sent  standard  deviation. 
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Anode  Current  Collector 


Fig.  7.  Schematic  of  electrolyte  filling  process  in  a  unit  cell. 

On  the  other  hand,  graphite  which  is  the  most  popular  anode 
material  undergoes  plastic  deformation  during  the  roll-pressing 
process  as  shown  in  Fig.  3.  Scanning  electron  microscopy  (SEM) 
images  of  the  cathode  and  anode  were  taken  after  compression, 
and  are  shown  in  Fig.  3.  SEM  examination  shows  that  the  cathode 
particles  maintained  a  spherical  shape,  whereas  the  anode  particles 
deformed  to  an  elliptical  shape  which  corresponds  to  the  results  of 
Wang  et  al.  [15  .  Therefore  the  deformation  of  anode  particle 
should  be  taken  into  consideration  for  the  simulation  of  wettability 
in  a  compressed  anode. 

For  the  reconstruction  of  porous  electrodes,  it  is  assumed  that 
the  electrode  consists  of  randomly  packed  spherical  particles  at 
arbitrary  positions.  The  structure  of  the  porous  electrode  is  con¬ 
structed  with  impermeable  solid  particles  and  pore  areas.  Spherical 
particles  are  generated  with  a  prescribed  radius  and  porosity.  Since 
the  newly  generated  particle  can  take  its  position  close  enough  to  a 
pre-existing  neighborhood,  overlapping  between  the  particles  is 
allowed.  Therefore,  the  actual  porosity  in  the  numerical  domain  is 
determined  by  the  ratio  of  the  number  of  void  grid  lattices  to  the 
number  of  grid  lattices.  It  is  assumed  that  compression  force  does 
not  destroy  the  solid  particles,  but  merely  reduces  the  pore  space. 
Based  on  these  assumptions,  the  microstructure  of  both  the  cath¬ 
ode  and  anode  are  constructed  as  shown  in  Figs.  4  and  5.  Figs.  4(a) 
and  5(a)  show  the  uncompressed  electrodes,  which  are  used  as  a 
baseline.  Four  different  compression  ratios  (10%,  20%,  30%,  and  40%) 
were  taken  into  consideration  at  the  cathode  and  anode  electrodes. 
As  the  compression  progress,  the  volume  of  the  particle  is 
conserved,  whereas  that  of  the  void  space  is  reduced.  This  indicates 
a  decrease  in  the  porosity.  As  previously  mentioned,  the  particle 
shape  at  the  cathode  remains  spherical  because  cathode  materials 
are  strong  in  compression,  and  the  particle  shape  at  the  anode 
becomes  elliptical  due  to  plastic  deformation.  In  the  compressed 
anode,  spherical  particles  are  compressed  in  a  through-plane  di¬ 
rection  and  extended  in  an  in-plane  direction. 

Fig.  6  shows  the  porosity  changes  in  the  cathode  and  anode  for 
various  compression  ratios.  Because  the  initial  porosity  was  set  to 
0.4,  the  porosity  is  expected  to  be  zero  when  the  compression  ratio 
is  over  40%.  On  the  other  hand,  the  predicted  porosity  is  greater 
than  zero.  This  is  because  the  reconstruction  is  made  in  the  two- 
dimensional  domain  by  allowing  an  overlap  between  the  parti¬ 
cles.  For  the  sake  of  clarity,  corresponding  error  bars  are  shown  in 


Fig.  9.  Liquid  electrolyte  distribution  in  the  cathode  when  the  compression  ratios  are 
(a)  0%,  (b)  10%,  (c)  20%,  (d)  30%  and  (e)  40%. 

Fig.  6.  The  figure  shows  that  the  porosity  in  the  cathode  is  larger 
than  that  in  the  anode.  This  is  due  to  the  particle  shape,  which 
remained  spherical  in  the  cathode,  but  deformed  to  an  elliptical 
shape  in  the  anode.  From  a  geometric  point  of  view,  the  void  space 
between  the  spherical  particles  is  greater  than  the  void  space  be¬ 
tween  the  elliptic  particles. 

23.  Numerical  procedure 

The  electrolyte  filling  process  in  a  unit  cell,  i.e.,  cathode,  anode, 
separator,  cathode  current  collector  and  anode  current  collector,  is 
shown  in  Fig.  7.  Based  on  results  reported  in  Ref.  [7  ,  the  wetting 
mechanism  is  described  as  follows.  First,  the  electrolytes  are 
transported  from  the  inlet  at  the  beginning  of  the  process.  Second, 
the  inflowing  electrolytes  penetrate  through  the  porous  electrodes 
and  separator  due  to  capillary  pressure.  Third,  the  electrolytes 
spread  into  the  electrodes  from  the  separator  because  the  elec¬ 
trolyte  transportation  in  the  separator  is  faster  than  that  in  the 
electrodes.  To  simulate  liquid  electrolyte  transport  in  a  porous 
electrode,  a  lattice  Boltzmann  (LB)  model  was  developed  as  shown 
schematically  in  Fig.  8.  To  simplify  the  LB  simulation,  the  compu¬ 
tational  domain  consists  of  an  electrode  and  a  separator.  The  two- 
dimensional  domain  size  is  1000  x  105  lu2  including  a  porous 
electrode  of  1000  x  100  lu2  and  a  separator  of  1000  x  5  lu2  (lu: 
lattice  unit,  1  lu  means  1  pm  of  the  physical  length  in  this  study).  It 
is  assumed  that  the  electrode  is  comprised  of  impermeable  solid 
particles  and  pore  areas,  and  the  separator  is  comprised  of  highly 
porous  material.  The  boundary  conditions  proposed  by  Zou  and  He 
[16],  who  implemented  no-slip  bounce  back  condition  on  the  walls, 
are  applied.  The  inlet  condition  was  set  on  the  left  boundary  with 
an  inflowing  velocity  of  7.0  x  10~5  lu.  The  resulting  Reynolds  and 
capillary  numbers  are  0.005  and  8.33  x  10  6,  respectively.  A  fixed 
pressure  condition  was  applied  at  the  right  boundary  and  the  top 
and  bottom  boundaries  were  set  to  a  no-slip  condition.  The  bounce 
back  condition,  i.e.,  the  incoming  distribution  function  is  simply 
equal  to  the  corresponding  outgoing  distribution  function  with  an 
opposite  momentum,  was  imposed  at  the  top  and  bottom  walls.  A 
porosity  of  0.4  and  a  particle  radius  of  5-7  lu  were  used  as  initial 
values  for  the  simulation  of  the  uncompressed  electrode.  Since  the 
roughness  can  change  the  contact  angle  of  particle  surface,  the 
neutral  wetting  case  with  a  contact  angle  of  90°  was  considered  in 
the  present  study.  The  finer  time  step  and  grid  size  ensure  better 


Fig.  8.  Computational  domain  of  two-dimensional  LB  model  with  boundary  conditions. 
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convergence  and  accuracy,  but  they  incur  a  high  calculation  cost. 
The  appropriate  time  step  (At  =  1  lu)  and  grid  size  (Ax  =  1  lu)  were 
selected  by  considering  the  motion  of  the  liquid  electrolyte  and 
gaseous  air.  The  computational  time  for  this  calculation  was  about 
4.3  x  105  s  per  9  x  107  time  steps  on  12  nodes  of  parallel  processors 
using  the  Linux  operating  system  with  a  3.30  GHz  clock  speed.  The 
validation  for  the  present  LB  model  was  reported  in  our  previous 
paper  [8],  so  it  will  not  be  repeated  here. 

3.  Results  and  discussion 

3  A.  Wettability  in  the  cathode 

Fig.  9  shows  the  liquid  electrolyte  distribution  in  the  cathode. 
Five  different  compression  ratios  (0%,  10%,  20%,  30%  and  40%)  are 
considered  for  a  lattice  time  of  6.0  x  107  lu.  Since  the  electrolyte 
distribution  did  not  show  any  difference  after  this  time,  the 
calculation  was  assumed  to  have  converged.  To  relate  the  physical 
space  to  the  lattice  space,  a  time  step  of  At  =  6.67  x  10~10  is  chosen. 
In  this  figure,  the  blue  and  light  gray  represent  the  liquid  electrolyte 
and  gaseous  air,  respectively.  The  liquid  electrolyte  driven  by  con¬ 
stant  velocity  at  the  inlet  pushes  gaseous  air  toward  the  outlet.  The 
electrolyte  transport  is  due  to  the  interplay  among  three  forces:  the 
pressure  force  FP  =  -A pA  from  the  connecting  electrolytes,  the 
capillary  force  from  the  surface  tension  Fs  =  -and,  and  the  resis¬ 
tance  force  Fr  from  the  existing  air  in  the  pores.  As  a  result,  the  air 
moves  driven  by  the  flow  motion  of  the  electrolyte.  A  more  detailed 
wetting  mechanism  can  be  found  in  our  previous  paper  [8  .  The 
results  of  the  present  study  show  that  the  higher  compression  ratio 
decreases  the  electrolyte  distribution  caused  by  lower  porosity.  In 
Fig.  9(a)  which  shows  an  uncompressed  cathode,  the  electrolytes 
are  distributed  in  large  areas.  Two  main  streams  are  observed 
which  are  transported  by  the  insertion  of  the  flow  from  the  inlet 
and  the  integration  of  the  flow  paths  spread  from  the  separator. 
These  are  distributed  mainly  at  the  inlet  and  the  bottom  of  the 
cathode.  Gaseous  air  is  observed  in  the  middle  and  at  the  outlet. 
The  former  is  trapped  by  the  electrolyte  and  the  latter  obstructs  the 
flow  path.  Fig.  9(b)  shows  a  10%  compressed  cathode,  and  the  area 
of  liquid  electrolyte  is  decreased  by  a  reduction  in  the  pore  of  the 
infiltrated  flow  paths.  Most  flow  paths  that  are  spread  from  the 
separator  do  not  reach  the  top.  This  is  attributed  to  the  increase  in 
the  blockage  of  the  flow  paths  and  the  resistance  force.  Fig.  9(c)  and 
(d)  shows  that  the  area  of  liquid  electrolyte  is  reduced  significantly 


Fig.  10.  Transient  saturation  profile  of  liquid  electrolyte  in  the  cathode  with  various 
compression  ratios. 


and  the  electrolyte  distribution  is  changed  when  compared  to 
Fig.  9(a)  and  (b).  This  is  because  the  electrolytes  selected  another 
preferential  path,  which  is  caused  by  an  increased  resistance  force. 
In  Fig.  9(e)  which  shows  a  40%  compressed  cathode,  the  electrolytes 
are  observed  mainly  at  the  bottom  and  are  rarely  found  at  the  top. 
This  is  due  to  the  increasing  number  of  disconnected  pore  throats 
resulting  from  the  blockage  of  particles  and  smaller  pores.  This 
indicates  that  the  electrolytes  permeate  slightly  into  the  pores  in  a 
highly  compressed  cathode. 

Fig.  10  shows  the  variations  of  liquid  electrolyte  saturation  in  the 
cathode  along  the  lattice  time.  A  comparison  is  made  with  various 
compression  ratios.  The  electrolyte  saturation  is  measured  using 
the  ratio  of  liquid  electrolyte  area  to  the  pore  space  area.  As  the 
lattice  time  increases,  electrolyte  saturation  increases  rapidly  until 
a  certain  lattice  time,  and  then  the  gradient  decreases.  When  the 
cathode  is  uncompressed,  the  electrolyte  saturation  increases 
steadily  due  to  the  creation  of  penetration  paths  and  the  extension 
of  the  paths  into  the  pores.  The  uncompressed  cathode  presents  the 
highest  electrolyte  saturation.  When  the  compression  ratio  is  10%, 
20%  and  30%,  the  electrolyte  saturations  decrease  by  approximately 
30-50%  with  a  gradual  gradient.  This  indicates  that  the  electrolyte 
has  several  penetration  flow  paths,  and  thus  the  electrolyte  trans¬ 
port  proceeds  slowly  despite  the  increased  resistance  force.  The 
electrolyte  saturation  has  several  irregular  slopes  at  1.4  x  107  and 
2.2  x  107  lu  for  CR  =  10%;  0.7  x  107  and  3.4  x  107  lu  for  CR  =  20%; 
and  0.5  x  107  lu  and  4.3  x  107  lu  for  CR  =  30%.  This  is  because  the 
flow  found  a  new  path  and  squeezed  the  electrolyte  into  the  pores. 
When  the  compression  ratios  are  40%,  the  electrolyte  saturation  is 
reduced  significantly  and  converges  in  the  early  stage,  indicating 
that  there  are  limited  penetration  flow  paths  in  the  cathode. 

Fig.  11  shows  the  liquid  electrolyte  saturation  profile  along  the 
cathode  thickness.  The  electrolyte  saturation  is  compared  with 
various  compression  ratios.  This  figure  shows  that  the  electrolyte 
saturation  decreases  along  the  thickness.  The  electrolyte  saturation 
is  observed  to  be  high  at  the  bottom,  but  this  becomes  lower  along 
the  thickness.  The  uncompressed  cathode  has  a  comparatively 
linear  decrease  along  the  thickness,  whereas  the  compressed 
cathode  presents  a  sharp  drop  near  the  bottom.  Because  the  highly 
compressed  cathode  has  a  high  resistance  force,  it  is  reasonable  to 
expect  lower  wettability  with  an  increasing  compression  ratio.  The 
irregular  slope  is  caused  by  non-uniform  air  distribution  that  is 
trapped  inside  the  cathode.  The  electrolyte  saturation  of  an 


Fig.  11.  Liquid  electrolyte  saturation  profile  along  the  cathode  thickness. 
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Fig.  12.  Liquid  electrolyte  distribution  in  the  anode  when  the  compression  ratios  are 
(a)  0%,  (b)  10%,  (c)  20%,  (d)  30%  and  (e)  40%. 

uncompressed  cathode  is  approximately  30%  at  the  top.  On  the 
other  hand,  the  compressed  cathode  presents  low  electrolyte 
saturation  in  this  region  which  is  attributed  to  the  electrolyte 
penetration  from  the  inlet.  This  indicates  that  the  electrolyte  flow 
paths  from  the  separator  can  hardly  reach  the  top  in  a  highly 
compressed  cathode.  This  result  confirms  that  the  thick  cathode 
with  high  compression  ratio  can  lead  to  low  utilization  of  the 
cathode  due  to  a  low  wetting  capability. 

32.  Wettability  in  the  anode 

Fig.  12  shows  the  liquid  electrolyte  distribution  in  the  anode 
when  the  lattice  time  is  6.0  x  107  lu.  For  comparison  with  elec¬ 
trolyte  distributions  in  the  cathode,  the  same  compressing  condi¬ 
tions  were  applied.  Fig.  12(a)  presents  the  electrolyte  distribution  in 
an  uncompressed  anode  which  is  used  as  a  baseline.  Fig.  12(b),  (c) 
and  (d)  shows  the  electrolyte  distribution  for  compression  ratios  of 
10%,  20%  and  30%,  respectively.  Compared  to  the  plots  shown  in 
Fig.  9,  the  through-plane  permeation  of  the  liquid  electrolyte 
decreased  and  the  flow  paths  of  the  branches  changed  their  di¬ 
rection  due  to  particle  deformation.  Because  the  particles  are 
compressed  in  a  through-plane  direction  and  are  expanded  in  an 
in-plane  direction,  the  in-plane  pore  size  is  reduced.  Since  creeping 
flow  through  a  small  pore  is  dominantly  influenced  by  capillary 
pressure,  and  since  the  electrolyte  selects  a  relatively  wide  throat, 
the  infiltrated  flow  path  stretches  to  the  in-plane  direction  rather 


Fig.  13.  Transient  saturation  profile  of  liquid  electrolyte  in  the  anode  with  various 
compression  ratios. 


Table  1 

Predicted  result  of  liquid  electrolyte  saturation. 


CR 

0% 

10% 

20% 

30% 

40% 

Cathode 

57.8 

40 

37.3 

29.6 

23.9 

Anode 

57.8 

34.3 

26.4 

19.9 

13.4 

than  the  through-plane  direction.  These  phenomena  can  be  clearly 
seen  in  Fig.  12(e)  which  shows  the  electrolyte  distribution  when  the 
compression  ratio  is  40%.  Because  the  compression  ratio  is  signifi¬ 
cantly  high,  pore  sizes  are  substantially  reduced  and  many  pore 
paths  are  blocked  by  the  particles.  As  a  result,  most  electrolytes  are 
observed  near  the  bottom. 

Fig.  13  shows  variations  of  liquid  electrolyte  saturation  in  the 
anode  along  the  lattice  time.  As  the  compression  ratio  increases, 
the  electrolyte  saturation  decreases.  When  the  compression  ratio  is 
10%,  the  electrolyte  saturation  is  reduced  to  one-half  in  comparison 
with  the  uncompressed  anode.  This  is  attributed  to  a  decrease  in 
the  through-plane  flow  path.  When  the  compression  ratios  are  20% 
and  30%,  the  electrolyte  saturation  increases  rapidly  at  2.0  x  107 
and  3.2  x  107,  respectively.  This  is  due  to  the  squeezing  of  the 
electrolyte  that  found  a  new  path  in  the  pores,  as  previously 
described.  When  the  compression  ratio  is  40%,  the  electrolyte 
saturation  is  reduced  significantly  and  converges  in  the  early  stage. 
This  is  caused  by  a  high  compression  ratio  that  gives  limited  inva¬ 
sion  routes,  and  thus  the  flow  path  could  not  find  the  linked  throat 
to  penetrate  into  the  pores.  Comparing  the  cathode  and  anode,  the 
cathode  has  higher  electrolyte  saturation  than  the  anode  at  each 
compression  ratio,  as  shown  in  Table  1.  Considering  that  the 
porosity  of  the  anode  is  slightly  lower  than  that  of  the  cathode  as 
shown  in  Fig.  6,  it  can  be  concluded  that  the  effect  of  particle  shape 
is  significant. 

Fig.  14  shows  the  liquid  electrolyte  saturation  profile  along  the 
anode  thickness.  The  electrolyte  saturation  is  decreased  along  the 
anode  thickness.  The  electrolyte  saturation  of  compressed  anodes 
decreases  rapidly  near  the  bottom  of  electrode.  This  increases 
slightly  at  30  lu,  and  drops  thereafter.  The  increase  of  the  electrolyte 
saturation  in  this  region  is  due  to  electrolyte  intrusion  from  the 
inlet  that  is  not  spread  from  the  separator,  as  shown  in  Fig.  12.  For 
all  the  compressed  anodes,  the  electrolyte  paths  from  the  separator 
do  not  reach  the  top.  This  suggests  that  the  resistance  force  is  very 
high  in  a  through-plane  direction.  Compared  to  the  electrolyte 


Fig.  14.  Liquid  electrolyte  saturation  profile  along  the  anode  thickness. 
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saturation  profile  in  the  cathode  shown  in  Fig.  11,  the  electrolyte 
transport  in  an  in-plane  direction  is  increased,  whereas  that  in  a 
through-plane  direction  is  decreased. 

The  simulation  results  confirm  that  the  compression  ratio 
significantly  affects  electrolyte  penetration  into  a  porous  anode, 
and  suggest  that  the  particle  shape  is  important  for  the  elec¬ 
trolyte  penetration  in  the  anode.  This  further  suggests  that  a  thin 
anode  is  favorable  for  an  electrode  design  with  high  packing 
density. 

4.  Conclusion 

The  wettability  of  the  cathode  and  anode  has  been  investigated 
using  two-dimensional  LBM  simulation.  The  LB  model  was 
developed  to  simulate  the  effect  of  the  compression  ratio  on 
wettability.  The  reconstruction  of  spherical  particles  in  the  cath¬ 
ode  and  elliptical  particles  in  the  anode  was  established  by 
assuming  that  the  cathode  materials  are  strong  in  compression 
and  the  anode  materials  undergo  plastic  deformation.  This  model 
enables  us  to  visualize  the  electrolyte  penetrating  process  through 
micro-pores,  and  compare  the  wetting  capability  between  the 
cathode  and  anode.  It  should  be  noted  that  the  electrolyte  trans¬ 
port  through  the  front  and  back  was  not  considered  due  to  the 
limitations  of  two-dimensional  simulation  developed  in  this  work. 
Extension  of  the  model  to  three-dimensional  may  improve  the 
accuracy,  but  this  is  beyond  of  the  scope  of  this  study.  The  results 
show  that  the  effect  of  the  compression  ratio  on  wettability  is 
significant.  The  area  of  liquid  electrolyte  in  the  compressed  elec¬ 
trodes  is  reduced  compared  to  uncompressed  electrodes.  As  the 
compression  ratio  increases,  the  wetting  capability  is  decreased. 
Because  the  anode  particles  have  a  through-plane  compression 
and  in-plane  expansion,  the  through-plane  permeation  of  liquid 
electrolyte  decreases.  This  causes  a  lower  wetting  capability  in  the 
anode,  which  results  in  low  utilization  of  the  anode.  On  the  other 
hand,  the  cathode  particles  maintain  a  spherical  shape  as  the 
compression  progresses.  This  merely  decreases  the  porosity  by 
providing  an  environment  such  that  infiltrated  flow  paths  can 
easily  spread  from  the  separator  and  penetrate  in  a  through-plane 
direction.  This  suggests  that  the  effect  of  particle  shape  on 
wettability  is  significant. 

The  results  of  the  present  study  provide  a  mechanism  for  elec¬ 
trolyte  transport  in  a  compressed  electrode.  The  electrolyte  trans¬ 
ports  into  the  porous  electrode  via  two  major  paths:  flow  from  the 
inlet  and  spread  from  the  separator.  Both  paths  depend  on  the  pore 
size.  Spread  from  the  separator  strongly  depends  on  the  particle 
shape.  Considering  that  the  size  of  an  EV  or  HEV  battery  is  very 
large,  spread  from  the  separator  is  more  effective  than  flow  from 
the  inlet.  Therefore  careful  control  of  the  particle  shape  and  elec¬ 
trode  thickness  is  needed. 
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Nomenclature 

A  area,  m2 

Ca  capillary  number 

CR  compression  ratio 

Cs  speed  of  sound,  m  s-1 

d  pore  diameter,  m 

ez  lattice  velocity,  m  s-1 

f  distribution  function 

Fads  adhesion  force 

f'ext  external  force 

Pint  interparticle  force 

Fp  pressure  force 

Fr  resistance  force 

Fs  capillary  force  from  the  surface  tension 

G  Green’s  function 

Se  electrolyte  saturation 

A p  pressure  difference,  Pa 

t  time,  sec 

u  velocity,  m  s-1 

w  weight  coefficient 

x  position  coordinate,  m 

Greek  letters 

v  kinetic  viscosity,  Pa  s 

a  surface  tension,  N  nrT1 

p  density,  kg  nrT3 

t  relaxation  time 

function  of  density 

Subscripts  and  superscripts 
e  electrolyte 

eq  equilibrium 
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